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A microscopic calculation of the reaction cross-section for nucleon-nucleus scattering has been 
performed by explicitly coupHng the elastic channel to all particle-hole (p-h) excitation states in the 
target and to all one-nucleon pickup channels. The p-h states may be regarded as doorway states 
through which the flux flows to more complicated conflgurations, and subsequently to long-lived 
compound nucleus resonances. Target excitations for '*°''**'Ca, ^^^Ni, ^"Zr and ^^^Sm were described 
in a QRPA framework using a Skyrme functional. Reaction cross sections calculated in this approach 
were compared to predictions of a fitted optical potential and to experimental data, reaching very 
good agreement. Couplings between inelastic states were found to be negligible, while the couplings 
to pickup channels contribute significantly. For the first time observed reaction cross-sections are 
completely accounted for by explicit channel coupling, for incident energies between 10 and 40 MeV. 



PACS numbers: 24.10.Eq, 25.40.-h, 24.50.-|-g 

A quantitative description of nucleon-nucleus reactions 
is crucial for a broad variety of applications, including as- 
trophysics, nuclear energy, radiobiology and space science 
[TJ H] . A fully microscopic description of such reactions 
is quite complex and resource-consuming, as one needs 
to consider not only the desired outcome in an exit chan- 
nel, but also the interference and competition with all 
other possible outcomes. A successful account of elastic 
nucleon-nucleus scattering, for example, has to include 
the effects from the excitation of non-elastic degrees of 
freedom, such as collective and particle-hole (p-h) excita- 
tions, transfer reactions, etc. Formally, these non-elastic 
effects can be accounted for by the projection-operator 
approach of Feshbach 13] . The picture that emerges is one 
in which flux is removed from the elastic channel by cou- 
plings to the non-elastic degrees of freedom. An optical 
potential can therefore be defined [3l |4] as the effective 
interaction in a single-channel calculation that contains 
the effects of all the other processes that occur during 
collisions between nuclei. Optical potentials play a very 
important role in the description of nuclear reactions. 
They are extensively used to describe the interactions of 
projectile and target in the entrance channel, and the in- 
teraction of ejectile and residual nuclei after the reaction; 
they are crucial ingredients in direct-reaction as well as 
statistical (Hauser-Feshbach) calculations. 

Most widely used are phenomenological optical poten- 
tials fitted to reproduce experimental data sets. They 
have been extremely successful for many applications in- 
volving nuclei in the range of the fits. Unfortunately, such 
adjustable potentials make strong assumptions about lo- 



cality and momentum dependence that are probably not 
justified. In addition, for nuclei lying outside the range 
of the fits, such as the nuclei produced at rare-isotope 
facilities, in the r-process, and in advanced reactor ap- 
plications, this can lead to unquantifiable uncertainties. 
To achieve a better understanding of nuclear reactions 
and structure it is important to calculate optical poten- 
tials by first-principle methods. 

Within microscopic reaction theory, an optical poten- 
tial is comprised of two components. The first is a real 
bare potential, the diagonal potential within the elastic 
channel, which is generally obtained by folding the nu- 
cleon distributions of both nuclei with a nucleon-nucleon 
effective interaction. The second is a complex dynamic 
polarization potential which arises from couplings to in- 
elastic states. The resulting optical potential is composed 
of an imaginary potential and a real part usually slightly 
different from the bare potential. The former gives rise 
to absorption of fiux from the elastic channel to the other 
reaction channels, and is hence directly connected with 
observed reaction cross-sections. 

Several attempts have been made to generate optical 
potentials from microscopic approaches. Some have used 
the so-called nuclear matter approach [5], which pro- 
vides accurate results at nucleon energies > 50 MeV 
[5]. Recently, new methods based on self-energy the- 
ory have been implemented [7], and new calculations, 
which combine a nuclear matter approach and Hartree- 
Fock-Bogoliubov (HFB) mean field structure model, pro- 
vide encouraging results for neutron scattering below 15 
MeV [H] . Earlier attempts used the nuclear structure ap- 



2 

















It 


I'l' 


• * • . 




T 


p+ Ca 


• J. F. Turner etal., 1964 
■ R. F. Carlson etal., 1975 

♦ J. F. DIcello et al., 1970 








— 1 1 1 1 — 


— 1 1 1 1 1 1 




- 






• R. F. Carlson et al., 1994 
— 1 1 1 1 1 1 


4 






— 1 — 1 — 1 — 1 — 

m 


* ^-m-f s- 








p+ Nl 


• J. J. Menetetal., 1971 
■ J. F. Turner etal., 1964 

♦ T. Eliyakut-Roshko et al., 1995 















10 15 20 25 30 35 40 



E,a, (MeV) 

FIG. 1: (color online). Total reaction cross-section as a func- 
tion of the incident energy for p -I- ''"Ca, p + "^^Ca and p 
_l_ 58-|>^j The results are shown for couplings to the inelas- 
tic states lying below 30 MeV (dashed lines), to the inelastic 
and transfer channels (dash-dotted lines) and to the inelas- 
tic and transfer channels with non-orthogonality corrections 
(solid lines). The Koning-Delaroche [T3] optical model cal- 
culations are shown as short-dashed lines. The lines serve as 
guidance to the eye as calculations were performed only for 
Slab = 10, 20, 30 and 40 MeV. Data from Refs. [T4HI9] . 

proach, which is more suitable at energies below 50 MeV 
[S], and calculated second-order diagrams using particle- 
hole propagators in the RPA approximation [51-fT5] . How- 
ever, these were not able to fully explain observed absorp- 
tion: e. g., in Ref. [H], the couplings could account only 
for « 44% of the nucleon-nucleus absorption and, in Ref. 
|12j . only for w 71% including charge exchange. 

In this Letter, we report on the first major step to- 
wards achieving a complete microscopic calculation of the 
reaction cross sections for both neutron and proton in- 
duced reactions on a variety of medium-mass targets. A 
summary of our results for protons on ""^'^^Ca and ^^Ni is 
given in Figure [T] where our best predictions are the solid 
lines. The other curves will be described below, where 
we discuss our method in more detail. 

Our calculations aim to obtain the reaction cross sec- 
tions of all the open channels that can be reached in one 
step from elastic scattering. We then use a doorway ap- 
proximation, which takes the total flux leaving the elastic 
channel to all possible first-order channels to be indepen- 
dent of what happens afterwards: a nucleon later might 



escape as a free nucleon, the flux might equilibriate to 
compound-nuclear resonances, etc. 

To generate sets of excited states, we use RPA and 
QRPA structure models for finite nuclei, which start from 
HFB structure models based on energy-density function- 
als. For each excited state, we calculate the one-body 
transition density and corresponding transition potential 
by the methods of [IHl [H], which we use within large 
coupled-channels calculations. In addition to inelastic 
excitations, we also include couplings to pickup channels. 

To obtain the initially occupied proton and neutron 
levels in a nucleus, we use the Skryme energy-density 
functional SLy4 22, Table 1], a parametrization designed 
to describe systems with arbitrary neutron excess, from 
stable to neutron matter, by improving isotopic prop- 
erties, which overcomes deficiencies of other interactions 
away from the stability line. A HFB calculation gives the 
particle and hole levels of a given nucleus and fixes the 
p-h basis states for generating excited states within the 
framework of (Q)RPA, thus accounting for correlations 
caused by the residual interactions within the target. 

Our scattering effective nucleon-nucleon interaction is 
of Gaussian shape, with parameters matched to the vol- 
ume integral and r.m.s. radius of the M3Y interaction at 
40 MeV; it includes a knock-on exchange correction [23] . 
In momentum space, the central effective interaction is 
v^il) = C(7r/Ai|)3/2e-«'/(2MT)^^ ^ith = -24.1921 
MeV and fig — 0.7180 fm~^ for the isoscalar part of the 
interaction and Vq^ = 11.3221 MeV and = 0.7036 
fm^^ for the isovector component. We do not include any 
imaginary part in this effective interaction, as our aim 
is to include all non-elastic excitations explicitly in our 
model. We convolute v'^{q) with the transition densities 
to generate the configuration-space transition potentials. 
The bare potential in the elastic channel is the single- 
folded potential using the ground-state density from the 
HFB calculation. For simplicity, this potential was also 
used for all excited states. 

To explore the relative importance of the various con- 
tributions to the reaction cross section, we carried out a 
series of calculations: 

(1) Inelastic coupled-channels calculations will be 
shown for reactions involving protons and neutrons scat- 
tered by the nuclei ^°Ca, ^^^Ca, ^^Ni, ^^Zr and "''Sm, 
coupling the ground state to all levels with excitation 
energy {E*) lying below some limit, according to the 
QRPA model. The QRPA states above the particle emis- 
sion threshold are used to approximate exact scattering 
waves. Recent studies have shown that such wave func- 
tions contain large density distributions outside the nu- 
clear radius [24] . When used in reaction calculations they 
accurately represent the continuum [25'. Thus, processes 
containing one nucleon in the continuum (plus the inelas- 
tically scattered projectile) are included in our model. 

(2) Additional couplings between excited states were 
considered as predicted by the RPA model. These cou- 
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plings, however, were found to be negligible for scattering 
energies above 10 MeV, allowing us to disregard them in 
the subsequent calculations. 

(3) Consideration of a finite-range interaction in a HFB 
description of the target structure. For reactions of nu- 
cleons scattered by ^"Zr, the reaction cross section results 
using the QRPA model with the SLy4 force were found 
to be practically equivalent to the results found using 
RPA states and transitions with the Gogny DIS force 
[26j . This was observed despite the proton pairing gap 
of 1.2 MeV of ^"Zr. 

(4) Inclusion of one-nucleon pick-up channels: pick-up 
channels play an important role in nucleon-nucleus scat- 
tering [571 [2H]- Coupled reaction channels (CRC) cal- 
culations were performed, including all the channels for 
the formation of a deuteron, picking up the appropriate 
nucleon from occupied levels in the target. For trans- 
fers, we approximate the HFB target states by bound 
single-particle states in a Woods-Saxon potential, with 
the radii fitted to reproduce the volume radii and Fermi 
energy obtained by the HFB calculations. The volume 
diffuseness and spin-orbit parameters were taken from 
Koning-Delaroche optical potentials [13] at i?iab — 0, 
with spin-orbit radii adjusted by the same factor used to 
fit the volume part to HFB radii. To overcome numeri- 
cal limitations, we coupled explicitly only to the transfer 
channels, incorporating all inelastic couplings in the in- 
elastic optical potential already calculated in (1). 

CRC calculations require, in addition to the scatter- 
ing potentials in the incoming channel, a scattering po- 
tential between the deuteron and the remaining target. 
We adopted the Johnson-Soper [21] prescription as it in- 
cludes the effects of deuteron breakup in adiabatic (sud- 
den) approximation. In this prescription, the deuteron 
potential is the sum of the individual neutron and pro- 
ton potentials with the target. For the real parts we used 
the diagonal transition potentials of the corresponding 
nucleon-nucleus reaction and, for the imaginary parts, 
the sum of the imaginary parts of the Koning-Delaroche 
[13j optical potential for protons and neutrons. That is, 
fitted parameters are used in the imaginary part of the 
deuteron potential, while we leave for future work to cal- 
culate deuteron and nuclear potentials self-consistently. 

To assess the success of our large-scale coupled- 
channels approach, we compare the calculated reaction 
cross section to that obtained by one of the best available 
phenomenological optical potentials, henceforth referred 
to as a^^ dSj. 

We examined the convergence with respect to maxi- 
mum excitation energy, and found that convergence of 
the inelastic calculations requires coupling of all excited 
levels below the scattering energy (i.e. all open channels). 
This behavior is observed for each partial wave as well 
as at all energies, as is iUustrated for p + "^"Ca in Figs. 
[2]and[3j respectively. 

Although the reaction cross section increases with the 
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FIG. 2: (color online). Reaction cross-section as a function of 
the partial wave for the reaction p -I- ^''Ca at i5iab ~ 20 MeV. 
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FIG. 3: (color online). Total reaction cross-section as a func- 
tion of the incident energy for the reaction p -f ^''Ca, for the 
different inelastic calculations. The short-dashed line shows 
the results using the Koning-Delaroche [13] optical potential. 
The lines serve as guidance to the eye as calculations were 
performed only for Siab = 10, 20, 30 and 40 MeV. 

number of coupled states, to the limit where all open 
channels are coupled. Figure |3] shows that these inelas- 
tic couplings account only for a small fraction (« 23% 
at £;iab = 30 MeV) of a^^. However, after including 
couplings to the pickup channels through the CRC cal- 
culations, a large increment is found, approaching cr{^^ 
and the experimental data, as can be seen in Figure[T] An 
even better agreement can be obtained after we include 
the non-orthogonality terms [301 P- 226] in the CRC cal- 
culations, also shown in Figure [l] This correction arises 
because at small radii the deuteron bound state is not 
orthogonal to bound states occupied in the target. 

This work focuses on reaction cross sections, which test 
the modulus of the S-matrix elements. Additional in- 
sights can be gained, e. g., from elastic angular distribu- 
tions. Preliminary calculations of these give reasonable 
agreement with measured cross sections. 

In Figure |4] we present the reaction cross sections ob- 
tained for nucleons scattered by the nuclei ^°^^*Ca, ^®Ni, 
^"Zr and ^^^Sm at an incident energy of 30 MeV, as a 
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FIG. 4: (color online). Normalized total reaction cross-section 
at iSiab = 30 MeV as a function of the area of the targets 
40,48^g^^ 58jg-j^ 902j. ^^(1 "*Sm, with ro = 1.2 fm. The hnes 
have the same meanings as in Figure [l] and the open symbols 
are the calculations for the different target nuclei studied. 
Filled symbols are experimental data from Refs. [15l ITTHTQ] 



function of the area of the target. The absorption is 
shown relative to the reaction cross-section of a black 
sphere, which is approximately the geometrical area of 
the target. It can be seen again that, despite the impor- 
tant contribution of all inelastic couplings to the reaction 
cross-section, a large amount of absorption is due to the 
pickup channels and the corresponding non-orthogonality 
corrections. Explicitly considering such couplings en- 
abled us to account for practically all of the non-clastic 
cross-sections in the studied reactions. 

In summary, we have calculated the reaction cross- 
sections for nucleon induced reactions on nuclei "'"^^^Ca, 
^^Ni, ^°Zr and ^^^Sm, by explicitly calculating the cou- 
plings to all the doorway transfer and (Q)RPA inelas- 
tic channels. We found that inelastic convergence is 
achieved when all open channels are coupled. While in- 
elastic couplings account for an important part of the 
reaction cross section, most contributions come from cou- 
plings to the deuteron pickup channel, in which case the 
non-orthogonality terms are significant. We obtain re- 
action cross sections that are in good agreement with 
phenomenological optical model results and experimen- 
tal data. Such results, using the doorway approximation, 
are an important milestone. Future work on couplings 
between different types of non-elastic processes will cal- 
culate higher-order corrections. 

This work represents the first complete microscopic 
calculation that uses basic interactions between nucle- 



ons within the nuclei to predict reaction observables for 
incident energy as low as 10 MeV. Using state-of-the- 
art nuclear structure models coupled with large-scale re- 
action computations allowed the accurate prediction of 
measurable quantities. This will serve as basis for fu- 
ture fully-consistent ab initio developments for a range 
of nuclei including unstable species. 

This work was performed under the auspices of the 
U.S. Department of Energy by Lawrence Livermore Na- 
tional Laboratory under Contract DE-AC52-07NA27344, 
and under SciDAC Contract DE-FC02-07ER41457. 
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